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Neutrino oscillations are solidly established, with a hint of CP violation just emerging. Similarly, 
there are hints of lepton universality violation in 6 —^ s transitions at the level of 2.6a. By assuming 
that the unitary transformation between weak and mass charged leptons equals the leptonic mixing 
matrix measured in neutrino oscillation experiments, we predict several lepton flavor violating (LFV) 

B meson decays. We are led to the tantalizing possibility that some LFV branching ratios for B 
decays correlate with the leptonic CP phase <5 characterizing neutrino oscillations. Moreover, we 
also consider implications for £i —>■ ijikik decays. 

PACS numbers: 


Introduction 


The historical discovery of the Higgs boson [1, 2] would 
have completed our picture of particle physics, were it 
not for the solid evidence we now have that neutrino fla¬ 
vors interconvert [3]. Apart from neutrino oscillations 
and cosmology, no other signs of new physics (NP) have 
been established. However, some indirect signs might 
have been found by the LHCb collaboration. In 2013, 
they have published the results of the measurement of 
a variety of observables in & —>■ s transitions. In some 
cases, the experimental result was found to be in clear 
tension with the Standard Model (SM) prediction. These 
include angular observables [4-7] in — ?► K*fjT [8], 
as well as a sizable suppression of several branching ra¬ 
tios [9, 10]. Recently, the LHCb announced new results 
based on the complete LHC Run I dataset ]11]. The in¬ 
clusion of new data has confirmed the anomalies, which 
are currently at the ~ 4cr level. Furthermore, in 2014, 
the LHCb collaboration found an intriguing indication of 
lepton universality violation in the ratio ]12] 


BR{B ^ KpL+pL-) 
BRiB ^ Ke+e-) 


0.745t!^;]]^!] ± 0.036. (1) 


This experimental measurement, obtained in the low 
dilepton invariant mass regime, is 2.6tT away from the 
SM result = 1.0003 ± 0.0001 ]13]. Although the 
statistical significance of this discrepancy is not enough 
to claim a discovery, it is highly suggestive that several 
independent global fits [14-18] have shown that this hint 
can be explained by the same type of new physics con¬ 
tributions as the previous 6 —> s anomalies. 

The violation of lepton universality usually comes to¬ 
gether with the violation of lepton flavor. Based on sym¬ 
metry arguments, Glashow, Guadagnoli and Lane [19] 
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recently argued that the observation of universality vio- 
lation in the lepton flavor conserving (LFC) B —> KCl 
decays implies the existence of the lepton flavor violating 
(LFV) processes B —>• KlfiJ (with i ^ j). The idea 
of LFV in B meson decays has been further explored 
in [20-24]. 

Here we take a step further in this direction. Since 
we lack a theory a flavor, we can not make definite pre¬ 
dictions for the LFV rates using the LFC ones as in¬ 
put. Hence we make the simplest alternative assump¬ 
tion, namely, that the unitary transformation between 
weak and mass charged lepton states is given by the lep¬ 
tonic mixing matrix measured in neutrino oscillation ex¬ 
periments. Under this assumption we make numerical 
predictions for several LFV observables in the B system. 
We emphasize that this assumption is not completely ad 
hoc. It will actually be a prediction in models where the 
leptonic mixing arises from the charged lepton sector. 
We refer to [25] for a general discussion and an example 
model. 


General aspects of the b ^ s anomalies 

The effective hamiltonian describing 6 —> s transitions 
can be expressed as: 

A 2 

- ^ ^ X (2) 

^ i 

Here Gf is the Fermi constant, e the electric charge and 
V the CKM matrix. The Wilson coefficients Ci and G' 
encode the different (SM and NP) contributions to the 
effective operators Oi and O'. The analysis of the avail¬ 
able experimental data on 6 —> s transitions reveals that 
the effective operators relevant for the resolution of the 
& —>■ s anomalies are: 

09 = 0^^ = (sjM) (/27» , (3) 

010 = O^o^ = (sjM) (/27“75M) , (4) 

where Pl = |(1 — 75) is the left-chirality projector. Sev¬ 
eral independent global fits [15-18] find a significant ten- 
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sion between the SM results for the Wilson coefficients 
of these operators and the experimental data. This can 
be clearly alleviated in the presence of NP contributions. 
According to the global fit [18], the Cg^ coefficient is the 
key to improve the fits. More precisely, one finds a rea¬ 
sonable agreement with data when NP provides a nega¬ 
tive contribution to Og^, with Cg^’^^ ~ —30%xCg^’^^. 
Similar improvements are found when NP enters in the 
SU{2)l invariant direction Cg^’^^ = —with 

C-MM.NP _ _^2% X 

Predicting lepton flavor violation in B meson decays 

Here we raise the following question: can the leptonic 
mixing matrix provide the required lepton flavor struc¬ 
ture in Og and Oio? And if so, what are the predictions 
for lepton flavor violation in the B sector? As suggested 
by global fits, let us assume that the relevant NP opera¬ 
tor contains a left-handed leptonic current. In this case, 
this operator can be generally written in the mass basis 
as: 

> (5) 

where 

= ClhcPLS, ( 6 ) 

Ji. = (7) 

and A is the energy scale of the NP inducing this oper¬ 

ator. The i,j indices denote the lepton flavor combina¬ 
tion characterizing the operator in eq. (7). The 3x3 
matrices C‘^ and completely determine the relations 
among the Wilson coefficients for different flavor choices. 
On the other hand, in the interaction (gauge) basis, O 
takes the same form, but the quark and lepton currents 
are written in terms of gauge eigenstates d! and P as 

Ji = ( 8 ) 

Jl, = ( 9 ) 

We now focus on the leptons. By combining eqs. (7) 
and (9) one finds the relation between and C^, 

C^ = UlC^Ue, (10) 

where Ui is the unitary matrix which relates the left- 
handed charged lepton gauge and mass eigenstates as 
I' = Similarly, the left-handed neutrino gauge and 
mass eigenstates are connected by another matrix, [7^, as 
v' = VuV. The product of these two matrices determines 
the leptonic charged current weak interaction, 

/ice = -^[W-fVT^LZ^'+h.c.] 

= -^^[iP-^ViCPL^ + ti.c.] , (11) 


where K = ujUi, is the leptonic mixing matrix measured 
in neutrino oscillation experiments. If C/i, =1, the left- 
handed neutrino gauge and mass eigenstates are the same 
and all the mixing is in the left-handed charged leptons. 
In this case K = Uj and eq. (10) leads to 

. ( 12 ) 

We do not attempt to give any model prediction for 
. Instead, we will assume that it is diagonal but with 
non-universal entries. In that case one can determine 
the required which, after using eq. (12), leads 
to a matrix compatible with the observations in 
b s transitions. In particular, the resulting must 
have a strong hierarchy between the ee and /r/r entries, 
P'ee ^ Order to induce a sizable correction to 

B —)■ kP'> fjT and a negligible one to i? —>■ . 


“Deriving” from neutrino oscillations 

Barring tuning of the parameters, we find two generic 
matrices in the gauge basis that lead to valid 
matrices in the mass basis. Their forms define our two 
scenarios: 

• Scenario A: = diag(0, e, 1) 

• Scenario B: = diag(e, 0,1) 

Here e <C 1 is a small parameter (interestingly enough, 
note that Ref. [19] considered = diag(0, 0,1), which 
corresponds to any of our scenarios in the limit e = 0) ^. 
Using the standard parameterization for the leptonic 
mixing matrix AT, one finds that in order to suppress 


the contributions to the 

ee Wilson coefficients. 

e must be 

close to 

eA = 

tan^ 013 

in scenario A, 

(13) 


eb = 

sin 012 
tan^ 013 

in scenario B. 

(14) 

COS^ 012 


Taking 3cr ranges for the mixing angles from the latest 
global fit to neutrino oscillation data [26], one finds the 
ranges [—0.10,-0.05] for scenario A and [-0.05,-0.03] 
for scenario B, irrespective of the neutrino mass spec¬ 
trum; normal and inverted hierarchies giving basically 
the same results. Interestingly, Oig ^ 0 implies e ^ 0, 
indicating a suggestive connection between quarks and 
leptons. 

We can now obtain for both scenarios. Let us first 
consider case A. Assuming e = ca and taking the best-fit 
values from [26], we find: 
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/ 0 -0.023+ 0.117e*^ 0.026 + 0.102 \ 

-0.023 + 0.117 e-*'5 0.005 cos <5 + 0.532 -0.001 cos (5 + 0.005 i sin <5 + 0.509 , (15) 

V 0.026+ 0.102 e-*'^ -0.001 cos <5 - 0.005 z sin (5 + 0.509 0.394 - 0.005 cos J / 


where <5 is the Dirac leptonic CP violating phase. In the Regarding case B, assuming now e = eb and taking the 
CP conserving case = 0) this matrix simplifies to best-fit values for the mixing angles from [26] , one finds 

/ 0 0.094 0.128 \ 

= 0.094 0.537 0.508 . (16) 

\ 0.128 0.508 0.389 / 




0 

0.011+ 0.117e-*^ 
- 0.012 + 0.102 e -*'^ 


0.011+ 0.117e*^ 
0.548-0.003 cos 5 
0.003* sin (5 + 0.489 


- 0.012 + 0.102 \ 
-0.003* sin (5+ 0.489 
0.003 cos ,5+ 0.416 / 


(17) 
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FIG. 1: The branching ratio of the decay B ^ Kejj, versus 
the CP violating phase S in scenarios A and B. The bands are 
obtained by taking the leptonic mixing angles within their la 
range w.r.t. the best-fit value (solid line) [26]. 



In the CP conserving case {6 = 0) this matrix simplifies 
to 


/ 0 
= 0.128 
\ 0.090 


0.128 0.090 \ 
0.545 0.489 
0.489 0.419 / 


(18) 


Lepton flavor violation in the B system 

The matrix can be used to make definite predic¬ 
tions for ratios of branching ratios in i? —>■ K£f £j decays, 


BR{B ^ KifiJ) 


2 Pnp ^ij 
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Here BR(H ^ KifiJ) = BR(H ^ KiJiJ) + BR(H ^ 
Ki-iJ) and BR(H ^ = (4.29 ± 0.22) x lO"'^ 

is the LHCb result [9], measured using the 3 fb“^ 
dataset after LHC Run I in the complete range, where 
= M2, is the dimuon invariant mass. The factor 
Pnp is the NP fraction of the B —>■ amplitude, 

Pnp = A1np/A 4 Total [19]. Using the results of the 
global fit [18], which gives ~ —12% x 

Pnp is found to be pnp —0.136 Finally, the 
factor accounts for phase space and charged lepton mass 
effects. These introduce sizable corrections for final 
states including r leptons. Using the results of Ref. [27], 
we find ~ 1 and ^re = — 0.63. Finally, we note 

that the parameterization in terms of pnp is only exact 
in the limit of vanishing non-factorizable contributions. 
However, we have found that these corrections are 
negligible for the processes we are interested in. 


Comparing the matrices for our two scenarios, we 
find that they are of the same of order of magnitude and 
the most significant difference lies in the terms involving 
6. This is what will allow us to relate B decays to the 
leptonic CP phase. 


For (5 = 0, we obtain the following predictions for the 


^ The authors of [19] derive their value for pNP from the LHCb 
RfC measurement, obtaining pNP ~ —0.159. 
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FIG. 2: Same as fig. (1) for the branching ratio of the decay 
B —> Kct. 


B K LFV transitions in scenario A, 


the analogous operators for the purely leptonic LFV pro¬ 
cesses —)■ tjik^k are: 

^4^ = ^ iml^PLiu) . (24) 


Here we assume that the scale of the NP responsible for 
the vectorial LFV operators is the same as the one rel¬ 
evant for B meson decays, eq. (5), although in full gen¬ 
erality these could be unrelated. The flavor structure of 
Ou = is given by the product which, fol¬ 

lowing the same prescription as for the B meson decays, 
can be written as = (kC^kA (kC^kA . 

\ ^ ij ^ ' 'mn 

The 04 e operator induces several ii —>■ £j£k(k decay 
processes: (i) £~ £j£k^k^ ^^id (ii) £~ (-'j^k^k (with 

k j). Their branching ratios can be written as [28] 


BR(£, ^ £,44) 


5127r3r, 


A4 ’ 


(25) 


BRiB ^ e [4.3, 6.2] X 10"^°, (20) 

BR(H ^ Ke^T^) G [0.4,1.3] x 10"^ , (21) 

BR(R ^ G [0.8,1.6] X 10"® . (22) 

These have been derived using the LHCb central value 

and taking the leptonic mixing angles in the preferred la 
ranges found by the fit [26]. The main generic prediction 
from our setup is thus 

BR(R ^ > BR(R ^ Ke^^i'f),BR{B Ke^r^) 

(23) 

However, experimentally the decay B —)■ is 

the easiest to search for and reconstruct. Indeed, elec¬ 
tron and tau final states are, 0(20%) and 0(80%) re¬ 
spectively, worse to reconstruct. Moreover future RUN 
II data will probe the region of 0(10“^®) for this channel 
providing a test of our scenario. 


where k = 2/3 when there are two identical leptons in the 
final state, and k = 1/3 otherwise, and m£. and T^. are 
the mass and decay width of the ii lepton, respectively. 
The coefficient Mijk takes the form in case (i), 

CtkC^k in case (ii). 

One can now use the experimental limits on these LFV 
branching ratios to derive bounds on A. Processes in¬ 
volving C'/"g are strongly suppressed and thus they do 
not provide meaningful bounds. This is the case of 
fjT e~e~e^, t~ — >■ e“e“e“'' and t~ — >■ /r“e“e+. 

•In contrast, the combined LHCb+BaBar+Belle limit 
BR(t“ — >■ /i“/i“/i+) < 1.2 X 10“® [29] translates into 
A > 6.7 TeV (in both scenarios, A and B). The other r 
decay modes lead to slightly less stringent bounds. Fu¬ 
ture B factories are expected to improve on the search 
for T~ —> fj,~ fi~ fj,'^ , with sensitivies to branching ratios 
as low as ^ 10“® [30], allowing us to probe NP scales up 
to A - 12 TeV. 


Rare B decays and leptonic CP violation 

One can now consider a scenario with a non-zero value 
of the CP violating phase 6 characterizing neutrino oscil¬ 
lations. In this case, we are led to the fascinating possi¬ 
bility that the LFV branching ratios for B meson decays 
will depend upon 6. Our results can be found in figs. (1) 
and (2) corresponding to the decay modes B —>■ 
and B —>■ Kr^e^ respectively. This would suggest an 
alternative way of probing 5 by using LFV B meson de¬ 
cays. 

ii —>• ijikt-k decays 

The same strategy can be extended to other LFV ob¬ 
servables if induced mainly by vectorial operators, as in 
eqs. (3) and (4). Assuming the same leptonic currents. 


Conclusions and discussion 

In summary, we have suggested that the universality 
and flavor violating 6 —> s anomalies may be related to 
the pattern of neutrino oscillations. By assuming that the 
unitary transformation between weak and mass charged 
lepton eigenstates is given by the leptonic mixing matrix 
measured in neutrino oscillations we predict several lep¬ 
ton flavor violating B meson decay rates. This way we are 
led to the thrilling possibility that some of the rare LFV 
B decay branching ratios correlate with the leptonic CP 
phase 6 that characterizes neutrino oscillations. Other 
lepton flavor violating processes processes such as ii —>■ 
ijikik have been considered in a similar manner. Im¬ 
proved measurements at Belle should probe new physics 
scale at the level A ~ 12 TeV. Relevant scenarios involve 
additional neutral currents, such as schemes containing 
an extra Z' boson with lepton universality violation in 
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B decays [31-34], or possibly some realizations of the 
electroweak symmetry SU{3)c SU{3 )l (S> U{l)x [35- 
40[. Such schemes should be taken seriously should the 
observed hints in the B sector persist. Finally, we note 
that in this paper we have focussed on the case where 
the violation of lepton universality is caused by NP in 
B —)■ with negligible contributions to B —>■ Kee. 

The alternative hypothesis is also plausible, though it 
has a lower constraining power since the electron chan¬ 
nel is experimentally somewhat less constrained than the 
muonic one. 


Note added 

A few days ago, an update of [18[ was presented in [41[. 
While this would change slightly the value of pmp used in 
our analysis, our main point remains and the numerical 
results are also left essentially unchanged. 
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